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We evaluate gravitational lensing as a technique for the detection of extrasolar moons. Since 2004 gravitational microlensing has been 
successfully applied as a detection method for extrasolar planets. In principle, the method is sensitive to masses as low as an Earth 
mass or even a fraction of it. Hence it seems natural to investigate the microlensing effects of moons around extrasolar planets. We 
explore the simplest conceivable triple lens system, containing one star, one planet and one moon. From a microlensing point of view, 
this system can be modelled as a particular triple with hierarchical mass ratios very different from unity. Since the moon orbits the 
planet, the planet-moon separation will be small compared to the distance between planet and star. Such a configuration can lead to a 
complex interference of caustics. We present detectability and detection limits by comparing triple-lens light curves to best-fit binary 
light curves as caused by a double-lens system consisting of host star and planet - without moon. We simulate magnification patterns 
covering a range of mass and separation values using the inverse ray shooting technique. These patterns are processed by analysing a 
large number of light curves and fitting a binary case to each of them. A chi-squared criterion is used to quantify the detectability of 
the moon in a number of selected triple-lens scenarios. The results of our simulations indicate that it is feasible to discover extrasolar 
moons via gravitational microlensing through frequent and highly precise monitoring of anomalous Galactic microlensing events with 
dwarf source stars. 

Key words, gravitational lensing - planetary systems - methods: numerical - methods: statistical 



1. Introduction 

By now hundreds of extrasolar planets have been detected 1 . For 
all we know, none of the newly discovered extrasolar planets 
offers physical conditions permitting any form of life. But the 
search for planets potentially harbouring life and the search for 
indicators of habitability is ongoing. One of these indicators 
might be the presence of a large natural satellite - a moon - 
which stabilises the rotation axis of the planet and thereby the 
surface climate (Benn 2001). It has also been suggested that a 
large moon itself might be a good candidate for offering habit- 
able conditions (Scharf 2006). In the solar system, most planets 
harbour moons. In fact, the moons in the solar system outnum- 
ber the planets by more than an order of magnitude. No moon 
has yet been detected around an extrasolar planet. 

The majority of known exoplanets has been discovered 
through radial velocity measurements, with the first success- 
ful finding reported by Mayor & Queloz (1995). This method 
is not sensitive to satellites of those planets, because the stellar 
"Doppler wobble" is only affected by the orbital movement of 
the barycentre of a planet and its satellites, though higher-order 
effects could play a role eventually. Here we consider Galactic 
microlensing, which has led to the discovery of several relatively 
low-mass exoplanets since the first report of a successful de- 
tection by Bond et al. (2004), as a promising technique for the 
search for exomoons. 

As early as 1999, it has been suggested that extraso- 
lar moons might be detectable through transit observations 
(Sartoretti & Schneider 1999), either through direct observa- 
tion of lunar occultation or through transit timing variations, 
as the moon and planet rotate around their common barycen- 
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tre, causing time shifts of the transit ingress and egress (cf. also 
Holman & Murray (2005)). In their simulations of space-based 
gravitational microlensing Bennett & Rhie (2002) mention the 
possibility of discovering extrasolar moons similar to our own 
Moon. Later that year, Han & Han (2002) performed a detailed 
feasibility study whether microlensing offers the potential to 
discover an Earth-Moon analogue, but concluded that finite 
source effects would probably be too severe to allow detections. 
Williams & Knacke (2004) published the quite original sugges- 
tion to look for spectral signatures of Earth-sized moons in the 
absorption spectra of Jupiter-sized planets. Cabrera & Schneider 
(2007) proposed a sophisticated transit approach using "mu- 
tual event phenomena", i.e. photometric variation patterns due 
to different phases of occultation and light reflection of planet 
and satellite. Han (2008) undertook a new qualitative study of a 
number of triple-lens microlensing constellations finding "non- 
negligible" light curve signals to occur in the case of an Earth- 
mass moon orbiting a 10 Earth-mass planet, "when the planet- 
moon separation is similar to or greater than the Einstein ra- 
dius of the planet". Lewis et al. (2008) analysed pulsar time-of- 
arrival signals for lunar signatures. Kipping (2009a,b) refined 
and extended the transit timing models of exomoons to include 
transit duration variations, and Kipping et al. (2009) examined 
transit detectability of exomoons with Kepler-class photometry 
and concluded that in optimal cases moon detections down to 
0.2 MEarth should be possible. 

We cover here several new aspects concerning the microlens- 
ing search for exomoons, extending the work of Han (2008). 
First, the detectability of lunar light curve perturbations is de- 
termined with a statistical significance test that does not need 
to rely on human judgement. Second, all parameters of the two- 
dimensional three-body geometry, including the position angle 
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of the moon with respect to the planet-star axis, are varied. Third, 
an unbiased extraction of light curves from the selected scenarios 
enables a tentative prediction of the occurrence rate of detectable 
lunar light curve signals. A more detailed account of this study 
is available as Liebig (2009). 

The paper is structured as follows: In Section 2, we recall 
the fundamental equations of gravitational microlensing relevant 
to our work. Our method for quantifying the detection rates for 
extrasolar moons in selected lensing scenarios is presented in 
Section 3. In Section 4, we discuss the astrophysical implica- 
tions of the input parameters of the simulations. Our results are 
presented in Section 5, together with a first interpretation and a 
discussion of potential problems of our method. 

2. Basics of gravitational microlensing 

The deflection of light by massive bodies is a consequence of the 
theory of general relativity (Einstein 1916) and has been experi- 
mentally verified since 1919 (Dyson et al. 1920), see Paczyriski 
(1996) for an introduction to the field or Schneider et al. (2006) 
for a comprehensive review. 

The typical scale of angular separations in gravitational lens- 
ing is the Einstein radius 9e, the angular radius of the ring of for- 
mally infinite image magnification that appears when a source at 
a distance D$ , a lens of mass M at a distance Dl, and the ob- 
server are perfectly aligned: 
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where G denotes the gravitational constant, c the speed of light. 
Dls is the distance between lens plane and source plane; in the 
non-cosmological distance scale of our galaxy Dls - Ds - Dl 
holds true. With the Einstein radius also comes the characteristic 
time scale of transient gravitational lensing events, the Einstein 
time fg = Dl0e/v ± with the transverse velocity v± of the lens 
relative to the source. 

Here we focus on the triple lens case with host star (S), 
planet (P) and moon (M), see Figure 1 for illustration. The lens 
equation can be expressed using complex coordinates, where r\ 
shall denote the angular source position and £ the image po- 
sitions, cf. choice of notation in Witt (1990) and Gaudi et al. 
(1998). stands for the angular position of the lensing body 
i. qjj is the mass ratio between lenses i and j (g, ; - = j^). The lens 

equation gets the following form, if the primary lens, the host 
star S , has unit mass and is placed in the origin of the lens plane, 
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This is a mapping from the lens plane to the source plane, which 
maps the images of a source star to its actual position in the 
source plane. As pointed out in Rhie (1997), and explicitly cal- 
culated in Rhie (2002), the triple-lens equation is a tenth-order 
polynomial equation in 

In microlensing, the images cannot be resolved. Detectable 
is only the transient change in magnitude of the source star, when 
lens and source star are in relative motion to each other. The total 
magnification y. of the base flux is obtained as the inverse of the 
determinant of the Jacobian of the mapping equation (2), 
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Fig. 1. Geometry of our triple-lens scenario, not to scale. Five 
parameters have to be fixed: The mass ratios qps = Mpi ane t/Ms ta i- 
and qmp = MMoon /^Pianet> the angular separations in the lens 
plane Ops and Qmp and the position angle of the moon (p. These 
parameters define uniquely the relative projected positions of the 
three bodies. 

Gravitational lensing changes the apparent solid angle of a 
source, not the surface brightness. The magnification fi is the 
ratio of the total solid angle of the images and the apparent solid 
angle of the unlensed source. 

For a status report of the past, present and prospective fu- 
ture of planet searching via Galactic microlensing the recent 
white papers of 2008/2009 are a good source of reference 
(Beaulieu et al. 2008; Dominik et al. 2008; Bennett et al. 2009; 
Gaudi et al. 2009). 

3. Method 

The simplest gravitational lens system incorporating an extraso- 
lar moon is a triple-lens system consisting of the lensing star, 
a planet and a moon in orbit around that planet, as sketched 
in Figure 1 . Most likely, lunar effects will first show up as no- 
ticeable irregularities in light curves that have been initially ob- 
served and classified as light curves with planetary signatures. 
To measure the detectability of a given triple-lens system among 
binary lenses, we have to determine whether the provided triple- 
lens light curve differs statistically significantly from binary-lens 
light curves. 

To clarify the terminology: We investigate light curves that 
show a deviation from the single-lens case due to the presence 
of the planetary caustic, i.e. which would be modelled as a star- 
plus-planet system in a first approximation. We call lunar de- 
tectability the fraction of those light curves which display a 
significant deviation from a star-plus-planet lens model due to 
the presence of the moon. We measure the difference between 
a given triple-lens light curve (which is taken to be the "true" 
underlying light curve of the event) and its best-fit binary-lens 
counterpart. The best-fit binary-lens light curve is found by a 
least-square fit. We then employe-statistics to see whether the 
triple-lens light curve could be explained as a normal fluctuation 
within the error boundaries of the binary-lens light curve. If this 
is not the case, the moon is considered detectable. 

Detection and characterisation are two separate problems in 
the search for extrasolar planets or moons, though, and we do 
not make statements about the latter. When characterising an ob- 
served light curve with clear deviations from the binary model, 
it is still possible that ambiguous solutions - lunar and non-lunar 
- arise. This does not reflect on our results. We simply give the 
fraction of triple-lens light curves significantly deviating from 
binary-lens light curves, without exploring whether they can be 
uniquely characterised. 

A qualitative impression of the lunar influence on the caus- 
tic structure can be gained from the magnification patterns 
in Figure 2. Extracted example light curves are presented in 
Figure 3. Going beyond the qualitative picture, in this paper we 
quantify the detectability of an extrasolar moon in gravitational 
microlensing light curves in selected scenarios. 
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Fig. 2. Details of the analysed triple-lens magnification maps for an example mass/separation scenario withes = 10~ 3 , quip = 10~ 2 , 
dps = 1.3 6e and 9mp = 1.0 (the standard scenario as summarised in Table 1), showing the caustic configurations for twelve 
different lunar positions completing a full circular orbit in steps of 30°. The relative positions of star, planet and moon are sketched 
in the lower left of each panel (not to scale). The side lengths of the individual frames are 0. 1 9e- A darker shade of grey corresponds 
to a higher magnification. The straight black lines mark source trajectories, the resulting light curves are displayed in Figure 3. 
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Fig. 3. Sample of triple-lens light curves corresponding to the source trajectories depicted in Figure 2. The solid light curves were 
extracted with an assumed solar source size, /^source = Rq', the thin, dashed lines show the same light curves with source = 3Rq. 
The magnification scale is given in negative magnitude difference (-Amag), i.e. the unmagnified baseline flux of the source is 0. 
The time scale in units of the Einstein time t E , assuming a uniform relative motion of source and lens. 
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3.1. Ray shooting 

The triple-lens equation (2) is analytically solvable, but 
Han & Han (2002) pointed out that numerical noise in the poly- 
nomial coefficients caused by limited computer precision was 
too high (~ 10~ 15 ) when solving the polynomial numerically 
for the very small mass ratios of moon and star (~ 10~ 5 ). To 
avoid this, we employ the inverse ray-shooting technique, which 
has the further advantage of being able to account for finite 
source sizes and non-uniform source brightness profiles more 
easily. It also gives us the option of incorporating additional 
lenses (further planets or moons) without increasing the com- 
plexity of the calculations. This technique was developed by 
Schneider & WeiG (1986), Kayser et al. (1986) and Wambsganss 
(1990, 1999) and already applied to planetary microlensing in 
Wambsganss (1997). 

Inverse ray-shooting means that light rays are traced from 
the observer back to the source plane. This is equivalent to trac- 
ing rays from the background source star to the observer plane. 
The influence of all masses in the lens plane on the light path is 
calculated. In the thin lens approximation, the deflection angle is 
just the sum of the deflection angles of every single lens. After 
deflection, all light rays are collected in pixel bins of the source 
plane. Thus a magnification pattern is produced (e.g. Figure 4). 
Lengths in this map can be translated to angular separations or, 
assuming a constant relative velocity between source and lens, 
to time intervals. The number of collected light rays per pixel is 
proportional to the magnification // of a background source with 
pixel size at the respective position in the source plane. The res- 
olution in magnification of these numerically produced patterns 
is finite and depends on the total number of rays shot. 

3.2. Light curve simulation 

We obtain simulated microlensing light curves with potential 
traces of a moon by producing magnification maps (Figure 4) of 
triple-lens scenarios with mass ratios very different from unity. A 
light curve is then obtained as a one-dimensional cut through the 
magnification pattern, convolved with the luminosity profile of a 
star. Only a few more assumptions are necessary to simulate re- 
alistic, in principle observable, light curves: angular source size, 
relative motion of lens and source, and lens mass. As a first ap- 
proximation to the surface brightness profiles of stars we use a 
profile with radius /^source and constant surface brightness (ig- 
noring limb-darkening effects). The light curve is sampled at 
equidistant intervals. In Section 4.7 we discuss physical impli- 
cations of the sampling frequency. In order to be able to make 
robust statistical statements about the detectability of a moon in 
a chosen setting, i.e. a certain magnification pattern, we anal- 
yse a grid of source trajectories that delivers about two hundred 
light curves, see also Figure 5. To have an unbiased sample, the 
grid is chosen independently of lunar caustic features, but all tra- 
jectories are required to show pronounced deviations from the 
single-lens light curve due to the planet. 

3.3. Fitting 

For each triple-lens scenario, we produce an additional magni- 
fication pattern of a corresponding binary-lens system, where 
the mass of the moon is added to the planetary mass. As a first 
approximation, one can compare two light curves with identi- 
cal source track parameters, see Figure 6(a). Numerically big 
differences can occur without a significant topological differ- 
ence, because the additional third body not only introduces ad- 




Fig. 4. Magnification pattern displaying the planetary caustic of 
a triple-lens system with mass ratios qps = 10~ 3 and q M p = 
10~ 2 , corresponding, e.g., to an M-dwarf, a Saturn-mass planet 
and an Earth-mass satellite of that planet. Darker shades corre- 
spond to areas of higher magnification. The third body shows 
its influence as the small perturbation on the bottom cusp of the 
planetary caustic. The position angle of the moon is <p — 90°, 
i.e. in the orientation of this map the moon is located above the 
planet. The complete set of parameters is given in Table 1. The 
side length of the pattern is 0.1 9e in stellar Einstein radii. A 
source with radius Rq (2Rq, 3Rq, 5Rq) is indicated as a black 
disc (dashed circles) in the lower left. Taking the straight line 
that cuts the lunar perturbation of the caustic as the source tra- 
jectory of a solar sized source results in the light curve shown 
in Figure 6 (solid line). Depending on the assumed source size, 
different light curves are obtained, see also Figure 9(a) to (c). 



ditional caustic features, it also slightly moves the location of 
the double-lens caustics. Since the perturbation introduced by 
the moon is small, however, we have good starting conditions 
using the light curve from the binary magnification pattern that 
corresponds to the triple-lens magnification pattern. Keeping the 
planet-star separation fixed and changing the planetary mass to 
the sum of the masses of planet and moon, we are likely to be 
close to the global best-fit values. In our simplified approach, 
we only vary the source trajectory to search for the best-fitting 
binary-lens light curve, which already yields satisfactory results, 
as can be seen in Figure 6(b). For the fitting procedure we use 
the simple, straight-forward and robust least-square method. 

3.4. Light curve comparison 

We want to use the properties of the x 1 -distribution for a test of 
significance of deviation between the two simulated light curves 
of the triple lens and the corresponding binary lens. 

Most commonly, the ^-distribution is used to test the 
goodness-of-fit of a model to experimental data. To compare two 
simulated curves instead, as in our case, one possible approach is 
to randomly generate artificial data around one of them and then 
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Fig. 5. Illustration of how the magnification pattern of a given 
star-planet-moon configuration is evaluated. In order to have an 
unbiased statistical sample, the source trajectories are chosen in- 
dependently of the lunar caustic features, though all light curves 
are required to pass through or close to the planetary caustic. 
This is realised through generating a grid of source trajectories 
that is only oriented at the planetary caustic. In the actual evalua- 
tion a 10 times denser grid than shown here was used. The mag- 
nification map has to be substantially larger than the caustic to 
ensure an adequate baseline for the light curve fitting. The side 
length of this pattern is 0.3 8 E , the position angle of the moon 
4> = 90°, all other parameters are chosen as listed in Table 1. 
Source sizes of 1 Rq, 2Rq, 3 Rq, and 5 Rq are indicated in the 
lower left. 



consecutively fit the two theoretical curves to the artificial data 
with some free parameters. Two '^-values' Q 2 and Q 2 ^ will re- 
sult. Their difference, AQ 2 = Q^-Q 2 , can then be calculated and 
a threshold value Ag t 2 hresh is chosen. The condition for reported 
detectability of the deviation is then AQ 2 > Ag 2 hresh . A detailed 
description of the algorithm and its application to microlensing 
was presented by Gaudi & Sackett (2000). 

We decided to use a different approach here: In particular, 
we were led to consider other possible techniques by the wish to 
avoid "feeding" our knowledge of the data distribution to a ran- 
dom number generator in order to get an unbiased and random 
^-distributed sample of Q 2 , when at the same time we have all 
the necessary information to calculate a much more representa- 
tive ^ 2 -value. We have developed a method to quantify signifi- 
cance of deviations between two theoretical functions that avoids 
the steps of data simulation and subsequent fitting. 

Put simply, we calculate the mean {Q 2 } of all '^-values' 
that would result, if data with a Gaussian scatter drawn from 
a triple-lens light curve is fitted with a binary light curve. If 
the triple-lens case and the best-fit binary-lens case are indistin- 
guishable, then Q 2 will be ^-distributed and (Q 2 ) will lie at or 
very close to (x 2 ) = number of degrees of freedom, the mean of 
the x 2 probability density function. If the triple-lens light curve 
and the binary-lens light curve differ significantly, then Q 2 is not 
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(a) Identical source trajectory. 
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(b) Best-fit source trajectory. 

Fig. 6. Triple-lens light curve (solid line) extracted from the 
magnification pattern in Figure 4, compared with a light curve 
extracted from the corresponding binary-lens magnification map 
(dashed line), with no third body and the planet mass increased 
by the previous moon mass. The difference of the two light 
curves is plotted as the dotted line. The central peak of the triple- 
lens light curve is caused by lunar caustic perturbations, that 
cannot easily be reproduced with a binary lens. In (a), identical 
source trajectory parameters are used to extract the binary-lens 
light curve and large residuals remain that can be avoided with 
a different source trajectory. In (b), the best-fit binary-lens light 
curve is shown (dashed) in comparison with the triple-lens light 
curve (solid). Here, the residuals can be attributed to the moon. If 
the deviation is significant, the moon is detectable in the triple- 
lens light curve. A source with solar radius at a distance of 8 kpc 
is assumed. 



X -distributed and (Q ) in particular will lie outside the expected 
range for ^-distributed random variables. In the latter case, the 
moon is detectable. This approach is presented in mathematical 
detail in Appendix A. 



4. Choice of scenarios 

This section presents the assumptions that we use for our sim- 
ulations. We discuss the astrophysical parameter space that is 
available for simulations of a microlensing system consisting of 
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star, planet and moon. By choosing the most probable or most 
pragmatic value for each of the parameters, we create a standard 
scenario, summarised in Table 1 and visualised in Figure 10, that 
all other parameters are compared against during the analysis. 

The observational search for microlensing events caused by 
extrasolar planets is carried out towards the Galactic bulge, 
currently limited to the fields of the wide-field surveys OGLE 
(Udalski et al. 1992) and MOA (Muraki et al. 1999). This leads 
to some natural assumptions for the involved quantities. The 
source stars are typically close to the centre of the Galactic bulge 
at a distance of D s = 8 kpc, where the surface density of stars 
is very high. For the distance to the lens plane we adopt a value 
of Di = 6 kpc. We assume our primary lens mass to be an M- 
dwarf star with a mass of Ms tar = 03Mq, because that is the 
most abundant type, cf. Figure 5 of Dominik et al. (2008). Direct 
lens mass determination is only possible if additional observ- 
ables, such as parallax, can be measured, cf. Gould (2009). We 
derive the corresponding Einstein radius of 9e(Ds , Di, Ms tar ) = 
0.32 milliarcseconds. 

Five parameters describe our lens configuration, cf. Figure 1: 
They are the mass ratios between planet and star q PS and be- 
tween moon and planet qMP, the angular separations 9p$ be- 
tween planet and star and 6 M p between moon and planet, and as 
the last parameter (f>, the position angle of the moon with respect 
to the planet-star axis. These parameters are barely constrained 
by the physics of a three-body system, even if we do require 
mass ratios very different from unity and separations that allow 
for stable orbits. 

The apparent source size Source and the sampling rate 
/sampled affect the shape of the simulated light curves. We have 
to define an observational uncertainty cr for the significance test 
in the final analysis of the light curves. 

4.1. Mass ratio of planet and star q PS 

We are interested in planets (and not binary stars). Accordingly, 
we want a small value for the mass ratio of planet and star 
Qps = ■ Our standard value will be 10~ 3 which is the mass 

'"Star 

ratio of Jupiter and Sun, or a Saturn-mass planet around an M- 
dwarf of 0.3 Mq. At a given projected separation between star 
and planet, this mass ratio determines the size of the planetary 
caustic, with a larger q PS leading to a larger caustic, qps is var- 
ied to also examine scenarios with mass ratios corresponding to 
a Jovian mass around an M-dwarf and to a Saturn mass around 
the Sun. Summarised: For the mass ratio between planet and star 
we use the three values q PS = 3.3 x 10~ 3 , 10 -3 , 3 x 10~ 4 . 

4.2. Mass ratio of moon and planet q M p 

To be classified as a moon, the tertiary body must have a mass 
considerably smaller than the secondary. The standard case in 
our examination corresponds to the Moon/Earth mass ratio of 
= 10~ 2 . We are generous towards the higher mass end, 

and include a mass ratio of 10 _1 in our analysis, corresponding 
to the Charon/Pluto system. Both examples are singular in the 
solar system, but we argue that a more massive moon is more 
interesting, since it can effectively stabilise the obliquity of a 
planet, which is thought to be favouring the habitability of the 
planet (Benn 2001). At the low mass end of our analysis we ex- 
amine q M p = 10~ 3 . In Figure 7, three different caustic interfer- 
ences resulting from the three adopted mass ratios are shown. 
Summarised: For the mass ratio between moon and planet we 
use the three values q M p = 10~ 3 , 10~ 2 , 10 -1 . 




Fig. 7. Caustic topology with different lunar mass ratios. From 
left to right, the lunar mass ratio increases from q M p = 10~ 3 to 
10 _1 . The side length of the displayed patterns is 0.05 9e- All 
other parameters as given in Table 1, in particular qps = 10~ 3 
and 8 M p = 1.0 0^. Source sizes of 1 Rq to 5Rq are indicated in 
the lower left. 



4.3. Angular separation of planet and star 6 PS 

The physical separation of a planet and its host star, i.e. the 
semi-major axis a PS , is not directly measurable in gravitational 
lensing. Only the angular separation 9 P s can directly be inferred 
from an observed light curve, where < 6p$Di < cips is valid 
for a given physical separation cips at a lens distance Di. The 
angular separation 8 PS of a binary will evoke a certain topology 
of caustics, illustrated e.g. in Figure 1 of Cassan (2008). They 
gradually evolve from the close separation case with two small 
triangular caustics on the far side of the star and a small central 
caustic at the star position, to a large central caustic for the inter- 
mediate case, when the planet is situated near the stellar Einstein 
ring, Ops =1.0 6e- This is also called the "resonant lensing" case, 
as described in Wambsganss (1997). If the planet is moved fur- 
ther out, one obtains a small central caustic and a larger isolated, 
roughly kite-shaped planetary caustic that is elongated towards 
the primary mass in the beginning, but becomes more and more 
symmetric and diamond-shaped if the planet is placed further 
outwards. 

In analogy, a "lunar resonant lensing zone" can be defined, 
when measuring the angular separation of moon and planet in 
units of planetary Einstein radii Q F E = -\fqp~s &e- Viewing our so- 
lar system as a gravitational lens system at Di — 6 kpc and with 
Ds = 8 kpc, Jupiter's Einstein radius in physical lengths would 
be 0.11 AU in the lens plane, that is almost 10 times the semi- 
major axis of Callisto. However, the smaller moons of Jupiter 
cover the region out to 2 Jovian Einstein radii. If there is a moon 
at an angular separation of the planet of 0.6 8 E < 9 M p ^ 3.0 
it will show its influence in any of the caustics topologies. But 
we focus our study on the wide-separation case for the follow- 
ing reasons: Regarding the close-separation triangular caustics, 
we argue that the probability to cross one of them is vanish- 
ingly small and, furthermore, the magnification substantially de- 
creases as they move outwards from the star. The intermediate 
or resonant caustic is well "visible" because it is always located 
close to the peak of the single-lens curve. But all massive bodies 
of a given planetary system affect the central caustic with mi- 
nor or major perturbations and deformations (Gaudi et al. 1998). 
There is no reason to expect that extrasolar systems are gener- 
ally less densely populated by planets or moons than the solar 
system. Since we are looking for very small deviations, identi- 
fying the signature of the moon among multiple features caused 
by low-mass planets and possibly more moons would be increas- 
ingly complex. The wide-separation case is most favourable, be- 
cause the planetary caustic is typically caused by a single planet 
and interactions due to the close presence of other planets are 
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highly unlikely (Bozza 1999, 2000). Any observed interferences 
can therefore be attributed to satellites of the planet. Planets with 
only one dominant moon are at least common in the solar sys- 
tem (Saturn & Titan, Neptune & Triton, Earth & Moon). Such a 
system would be the most straightforward to detect in real data. 

This makes a strong case for concentrating on the wide- 
separation planetary caustic. As the standard we choose a sepa- 
ration of Ops = 1.3 Be and also test a scenario with 9 PS = 1 .4 9 E . 
For comparison: the maximum projected separation of Jupiter at 
the chosen distances D$ and D E corresponds to 1 .5 solar Einstein 
radii. Summarised: For the angular separation between planet 
and star we use the two values PS = 1.36e, \A9 E . 

4.4. Angular separation of moon and planet 9 MP 

As the moon by definition orbits the planet, there is an upper 
limit to the semi-major axis of the moon omp- It must not exceed 
the distance between the planet and its inner Lagrange point. 
This distance is called Hill radius and is calculated as 

/ M P y 



for circular orbits, with the semi-major axis cips of the plane- 
tary orbit. This can be translated to "lensing parameters", but the 
equation changes into the inequality 

'"Hill > 6 PS D L ^qp S ^ , (4) 

because for a physical distance a PS , the apparent (i.e. projected) 
separation lies in the range < 9psD E < ap$ depending on the 
inclination towards the line of sight. Thus, the Hill radius con- 
straint a M p < '"Hill cannot be regarded as a strict limit. Though, 
turning this argument around, a minimum region of secured sta- 
bility exists for given mass ratios and angular planet-star separa- 
tion Ops . Bodies in prograde motion with an orbit below 0.5 rem 
can have long term stability, for retrograde motion the limit is 
somewhat higher at 0.75 r H iii (see Domingos et al. (2006) and 
references therein, particularly Hunter (1967)). The region of se- 
cured stability for our standard scenario has a radius of 

9 MP < 0.5 x 9 PS D L l^qpsj = 0.045 9 E = 1.53 9 P E , 

which leads to a M p = 0.09AU in our setting. Our choice of 
9 mp is also motivated by a desire to have caustic interactions 
between lunar and planetary caustic, a resonance, because only 
under this condition the moon will give rise to light curve fea- 
tures that are distinctly different from those of a low-mass planet. 
If the moon is located at an angular separation from the planet 
of - depending on the position angle of the moon - roughly 
0.6 9 F E < 9 mp 55 3.0 9 E , it will influence the planetary caus- 
tic. For angular separations larger than 3.0 9 E , the lunar caustic 
barely influences the planetary caustic, and if it is detected, it 
will be detected as a distinct secondary caustic. Han (2008) al- 
ready pointed out that the lunar signal does not vanish with an 
increased separation from the planet. But the lunar mass might 
be wrongly identified as a second, lower-mass planet rather than 
a satellite of the first planet. For comparison, the major moons 
of the solar system are all at considerably smaller angular sep- 
arations than 1 .0 9 F E , if observed as a gravitational lens system 
from D E — 6 kpc with Ds — 8 kpc, with our Moon being most 
favourable at 0.4 of Earth's Einstein radii at maximum projected 




Fig. 8. Magnification patterns with side lengths of 0.05 9 E show- 
ing the lunar perturbation of the planetary caustic for in- 
creasing angular separations of moon and planet of 9mp = 
0.8 9 E , 1 .0 9 F , 1 .2 9 E , while all other parameters are as in Table 1 . 
The topology of interaction varies with the angular separation. 
Source sizes of 1 Rq to 5 Rq are indicated in the lower left. 



separation. But as we have learnt through the discovery of hun- 
dreds of exoplanets, we do not have to expect solar system prop- 
erties to be mirrored in exoplanetary systems. 

We focus on configurations with strong interferences of plan- 
etary and lunar caustic, what one might call the "resonant case 
of lunar lensing", examples of which are displayed in Figure 8. 
Therefore, our standard case will have an angular separation of 
moon and planet of one planetary Einstein radius. In total, we 
evaluate four settings: 9 M p = 0.8#£, 1.06^, l.29 F , \A9 P E , in the 
chosen setting (Table 1), this corresponds to physical projected 
separations from 0.5 to 0.8 AU. 

4.5. Position angle of moon cf> 

The fifth physical parameter for determining the magnification 
patterns is the position angle of the moon <p with respect to 
planet-star axis, as depicted in Figure 1 . It is the only parame- 
ter that is not fixed in our standard scenario. Instead, we eval- 
uate this parameter in its entirety by varying it in steps of 30° 
to complete a full circular orbit of the moon around the planet. 
By doing this, we are aiming at getting complete coverage of a 
selected mass/separation scenario. Figure 2 visualises how the 
position of the moon affects the planetary caustic. It is not to be 
expected that we will ever have an exactly frontal view of a per- 
fectly circular orbit, but it serves well as a first approximation. 
Summarised: we choose 12 equally spaced values for the posi- 
tion angle <f> of the moon relative to the planet-star axis, and the 
results for the various parameter sets are averaged over these 12 
geometries. 

4.6. Source size Rsource 

The source size strongly influences the amplitudes of the light 
curve fluctuations and the "time resolution". Larger sources blur 
out finer caustic structures, compare Figure 9. For our source 
star assumptions, we have to take into account not only stel- 
lar abundances, but also the luminosity of a given stellar type. 
Main sequence dwarfs are abundant in the Galactic bulge, but 
they are faint, with apparent magnitudes of ~ 15 mag to 25 
mag at Ds = 8 kpc. In the crowded microlensing fields, they 
are difficult to observe with satisfying photometry by ground- 
based telescopes, even if they are lensed and magnified. Giant 
stars, with apparent baseline magnitudes ~ 13 mag to 17 mag, 
are more likely to allow precise photometric measurements from 
ground and, therefore, today's follow-up observations are biased 
towards these source stars. 



C. Liebig & J. Wambsganss: Exomoons as gravitational lenses 



9 



The finite source size constitutes a serious limitation to the 
discovery of extrasolar moons. In fact, Han & Han (2002) stated 
that detecting satellite signals in the lensing light curves will 
be close to impossible, because the signals are smeared out by 
the severe finite-source effect. They tested various source sizes 
(and planet-moon separations) for an Earth-mass planet with a 
Moon-mass satellite. They find that even for a KO-type source 
star (^source = 0.8 Rq) any light curve modifications caused by 
the moon are washed out. Han (2008) increased lens masses and 
by assuming a solar source size /^source = Rq, he finds that "non- 
negligible satellite signals occur" in the light curves of planets 
of 10 to 300 Earth masses, "when the planet-moon separation is 
similar to or greater than the Einstein radius of the planet" and 
the moon has the mass of Earth. We use a solar sized source 
as our standard case, and present results for four more source 
radii: Source = IRq, 2Rq, 3Rq, 5Rq, IORq. In terms of stellar 
Einstein radii this corresponds to five different source size set- 
tings from 1.8 x 10~ 3 9 E to 18 x 10~ 3 E . 

The ability to monitor smaller sources - dwarfs, low-mass 
main sequence stars - with good photometric accuracy is one 
of the advantages of space-based observations. Hence a satellite 
mission is the ideal tool for this aspect of microlensing, rou- 
tine detections of moons around planets can be expected with a 
space-based monitoring program on a dedicated satellite. 

4.7. Sampling rate /sampled 

Typical exposure times for the microlensing light curves are 10 
to 300 seconds at a 1.5m telescope. The frequency of observa- 
tions of an individual interesting event varies between once per 2 
hours to once every few minutes, for a single observing site only 
limited by exposure time and read-out time of the instrument. 
One example of this is the peak coverage of MOA-2007-BLG- 
400, see Dong et al. (2009), where frames were taken every six 
minutes. Higher observing frequency equals better coverage of 
the resulting light curve. A normal microlensing event is seen as 
a transient brightening that can last for a few days to weeks or 
even months. A planet will alter the light curve from fraction of 
a day to about a week; the effect of a low-mass planet or a moon 
will only last for a few hours or shorter, see e.g. Beaulieu et al. 
(2006). This duration is inversely proportional to the transverse 
velocity v ± with which the lens moves relative to the line of sight 
or the relative proper motion fi ± between the background source 
and the lensing system. We fix this at a typical v ± = 200 km/s 
for D L = 6 kpc or fi ± = 7.0 mas/year. 

For the simulations, a realistic, non-continuous observing 
rate is mimicked by evaluating the simulated light curve at in- 
tervals that correspond to typical frequencies, but allowing for 
a constant sampling rate, which in practise is often prohibited 
by observing conditions. We sample the triple-lens light curve 
in equally sized steps. With the assumed constant relative veloc- 
ity, this translates to equal spacing in time. Aiming for a rate of 
/sampled - is we choose a step size of 6 x 10~ 4 6 E . To see 

the effect of a decreased sampling rate, we have also examined 
the standard scenario (see Table 1) with sampling rates of factors 
of 1 .5 and 2 longer. Summarised: we use three different sampling 
rates, /sampled = 1/15 min, 1/22 min, 1/30 min. 



4.8. Assumed photometric uncertainty cr 

The standard error cr of a data point in a microlensing light 
curve can vary significantly. Values between 5 and 1 00 mmag 
seem realistic, see, e.g., data plots in the recent event analyses 
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Fig. 9. Illustration of the increasing source radius, compare 
Figure 6(b): (a) R SOUI ce = 2/?©, (b) 3 Rq, and (c) 5 Rq, all other 
parameters as given in Table 1 . The triple-lens light curve (solid) 
is obtained from the magnification pattern in Figure 4, the best- 
fit binary-lens light curve (dashed) from the corresponding mag- 
nification pattern without a moon. The difference between the 
two curves is plotted as a dotted line. 
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of Gaudi et al. (2008) or Janczak et al. (2009). The photometric 
uncertainty directly enters the statistical evaluation of each light 
curve as described in appendix A. We have drawn results for 
an (unrealistically) broad cr-range from 0.5 mmag to 500 mmag. 
An error as small as cr = 0.5 mmag for a V — 12.3 star has 
recently been reached with high-precision photometry of an ex- 
oplanetary transit event at the Danish 1 .54 m telescope at ESO 
La Silla (Southworth et al. 2009), which is one of the telescopes 
presently used for Galactic microlensing monitoring observa- 
tions, but this low uncertainty cannot be transferred to standard 
microlensing observations, where crowded star fields make e.g. 
the use of defocussing techniques impossible. In our discussion, 
we regard only the more realistic range from cr = 5 to 100 mmag. 
We adopt an ideal scenario with a fixed cr-value for all sampled 
points on the triple-lens light curve and set cr = 20 mmag to be 
our standard for the comparison of different mass and separa- 
tion scenarios. Summarised: we explore 5 values as photometric 
uncertainty: cr = 5, 10, 20, 50, 100 mmag. 



Parameter 


Standard value 


Comment 




D s 


8kpc 


distance to Galactic bulge 




D L 


6kpc 






Mstar 


0.3 Mq 


most abundant type of star 




V-x 


7 mas/year 


= v ± = 200 km/s at D L = 6 kpc 


* 


qps 


io- 3 


Jupiter/Sun mass ratio 


* 


Ops 


1.3 e E 


wide separation caustic 


* 


Imp 


io- 2 


Moon/Earth mass ratio 


* 


6mp 


i.o eg 


planetary Einstein radius 


* 


^Source 


R e 


brightness requirements vs. stellar 
abundance 


* 


./sampled 


1 frame 
15 minutes 


high-cadence observation 






20 mmag 


typical value in past observations 



Table 1. Parameter values of our standard scenario. Parameters 
marked with an asterisk (*) are varied in our simulations in order 
to evaluate their influence on the lunar detection rate and to com- 
pare different triple-lens scenarios. The fixed parameters lead to 
values for the Einstein ring radius, 9e = 0.32 mas, i.e. 1.9 AU 
in the lens plane, and the Einstein time, t E - 17 days. The 
lensed system is a Saturn-mass planet at a projected separation 
of 2.5 AU from its 0.3 Mq M-dwarf host, the Earth-mass satel- 
lite orbits the planet at 0.06 AU, i.e. 0.01 mas angular separation, 
cf. Figure 10. 



5. Results: detectability of extrasolar moons in 
microlensing light curves 

The numerical results of our work are presented and a first inter- 
pretation is made. We start by considering the result for a single 
magnification pattern in Section 5.1. We then present results that 
cover specific physical triple-lens scenarios (Section 5.2). We 
also evaluate the parameter dependence of the detectability rate, 
by varying each parameter separately. 

5.1. Evaluation of a single magnification pattern 

As described in Section 3, all magnification maps are evaluated 
by taking a large sample of unbiased 2 light curves and then de- 

2 I.e. independent of lunar features, but required to pass close to or 
through the planetary caustic. 



termining the fraction of light curves with significant lunar sig- 
nals. The results for a single magnification map (here the pattern 
from Figure 5), are displayed in Table 2. The percentages can 



cr in mmag 


Proportion of signifi- 
cantly deviating light 
curves 


Detectability 


5 


190/228 


83.3 % 


10 


132/228 


57.9 % 


20 


77/228 


33.8 % 


50 


39/228 


17.1 % 


100 


1/228 


0.4% 



Table 2. Detectability of the third mass in the magnification pat- 
tern of Figure 5, for different assumed photometric uncertainties. 



be interpreted as the "detection probability" of the moon in a 
random observed light curve displaying a signature of the plan- 
etary caustic - ignoring the possible complication to fully char- 
acterise the signal. In this map, the detectability of the moon is 
about one third, provided we have an observational uncertainty 
of cr = 20 mmag. Obviously, with a larger photometric error, one 
is less sensitive to small deviations caused by the moon. 

5.2. Results for selected scenarios 

In order to be able to make general statements about a given 
planetary/lunar system, we average the results for 12 mag- 
nification patterns representing a full lunar orbit with <f> = 
0°, 30°,..., 330°. 

5.2.1 . Variation of the photometric uncertainty 

The scenario described in Table 1 is evaluated for different pho- 
tometric uncertainties by analysing the 12 magnification maps 
shown in Figure 2. The results are displayed in Table 3. For 



cr in mmag Detectability 



5 


90.2 


% 


10 


59.8 


% 


20 


30.6 


% 


50 


12.8 


% 


100 


1.5 


% 



Table 3. Detectability of the moon in our standard scenario as a 
function of assumed photometric uncertainty. 



an photometric uncertainty of cr = 20mmag, the result is that 
about 30 % of the light curves show significant deviations from 
a best-fit binary-lens light curve and display detectable signs of 
the moon. We recall again our requirements: 

- The host planet of the satellite has been detected. 

- A small source size, up to a few solar radii, is required. 

- The moon is massive compared to satellites in the solar sys- 
tem. 

- We have assumed light curves with a constant sampling rate 
of one frame every 15 minutes for about 50 to 70 hours. 
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Fig. 10. Visualisation of the approximate lunar resonance zone (dark-grey ring), where planetary and lunar caustic overlap and 
interact. Also depicted is the planet lensing zone (light grey), which is the area of planet positions for which the planetary caustic(s) 
lie within the stellar Einstein radius. Distances and caustic inset are to scale. The physical scale depends not only on the masses 
of the lensing bodies, but also on the distances to the lens system and the background star. In our adopted standard scenario, cf. 
Table 1, the Einstein radius G E corresponds to 1.9 AU in the lens plane. This places the planet at a projected orbit radius of 2.5 AU, 
whereas the moon's orbit has a radius of 0.06 AU (corresponding to one planetary Einstein radius). The lunar resonance zone covers 
lunar orbit radii from 0.04 AU to 0. 1 8 AU. 



As would be expected, the detection rate increases with higher 
photometric sensitivity and decreases with a larger uncertainty. 
At a photometric error of 20 mmag or worse, the detection of a 
moon depends on whether the source trajectory passes through 
or close to the lunar caustic. The finite magnification resolution 
of the magnification patterns induces an uncertainty in the de- 
tectability results that increases with smaller cr-values, but it is 
negligible for the examined range. 

5.2.2. Changing the angular separation of planet and moon 

Table 4 lines up the changing detectability for different values of 
the projected planet-moon separation Omp, with all other param- 
eters as in standard scenario (Table 1). With the chosen range 



9 M p in 6 P E 


Detectability 


0.8 


29.6 % 


1.0 


30.6 % 


1.2 


30.4 % 


1.4 


28.1 % 



Table 4. Detectability of the moon as a function of the planet- 
moon separation for our standard scenario. 



from 6 mp = 0.8 Q P E to 1 .4 Q F E , we are within the regime of caustic 
interactions of planetary and lunar caustic. We find that the per- 
centage of detectability is almost constant at a 30 % level. Han 
(2008) defines regions of possible satellite detections as the re- 
gion of angular moon-planet separations between the lower limit 
of one planetary Einstein radius and an upper limit at the pro- 
jected Hill-radius average. From our results, we conclude that 
detections with 6 M p < 1 .0 ^ are well possible. We have not 
fully probed the lower limit region, but for 6mp < 0.5 Q P E the size 
of the lunar caustic perturbation decreases substantially. 

5.2.3. Variation of the moon mass 

Table 5 shows the changing detectability for a varying mass ratio 
qMP of lunar and planetary mass, with all other parameters as in 



our standard scenario (Table 1). Over the range of the moon- 



Imp 


Detectability 


io- 3 


2.4% 


io- 2 


30.6 % 


io- 1 


85.6 % 



Table 5. Detectability of the moon depending on moon-planet 
mass ratio. 



planet mass ratio, q M p = 10~ 3 , 10~ 2 , and 10~' the detection rate 
increases from 2.4 % to 85.6 %, compare also Figure 7. Not 
surprisingly: the more massive the moon, the easier it is to detect. 

5.2.4. Changing the planetary mass ratio 

Table 6 lists the changing detectability for a varying mass ratio 
between planet and star q PS . We adjusted the parameters of the 



qps 


Detectability 


3 x 10~ 4 


19.6 % 


1 X IO" 3 


30.6 % 


3 x IO" 3 


42.7 % 



Table 6. Dependence of lunar detectability as a function of the 
planet-star mass ratio, all other parameters as in our standard 
scenario. 



grid of source trajectories, so that the grid spacing scales with 
the total caustic size. For a decreasing mass of the planet, the ap- 
parent source size increases relative to the caustic size, so finer 
features will be blurred out more in the case of a smaller planet. 
Also the absolute change in magnification scales with the plan- 
etary mass. As expected, the detection efficiency is roughly pro- 
portional to the planetary caustic size. 
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5.2.5. Different sampling rates 

Table 7 presents the changing lunar detectability for different 
sampling rates. A lower sampling frequency lowers the detection 



Sampling step 


/sampled P^F 


Detectability 


size in 6 E 


minutes 




6 x 10- 4 


1/15 


30.6 % 


9 x 10- 4 


1/22 


26.8 % 


12 x 1(T 4 


1/30 


24.9 % 



Table 7. Detectability of the moon depending on the sampling 
rate of the simulated light curves. 



probability, but the effect is less pronounced than expected. This 
means it may be favourable to monitor a larger number of plan- 
etary microlensing events with high sampling frequency, rather 
than a very small number of them with ultra-high sampling. Our 
assumed sampling rates can easily be met by follow-up obser- 
vations, as they are presently performed for anomalous Galactic 
microlensing events. 

5.2.6. Source size variations 

We analysed simulated light curves of the standard scenario with 
five different source star radii. They cover the range between a 
star the size of our Sun and a small giant star with R Source = 
lOR Q atD s = 6kpc. 



^Source ln < 


h 


^Source in Rq 


Detectability 


1.8 X 10 


-3 


1.0 


30.6 % 


3.6 X 10 


-3 


2.0 


20.9 % 


5.4 x 10- 


-3 


3.0 


18.1 % 


9.0 x 10 


-3 


5.0 


4.1 % 


18.0 x 10 


-3 


10.0 


0.0% 



Table 8. Detectability of the moon depending on the size of the 
source star, all other parameters as in our standard scenario. 



A larger source blurs out all sharp caustic-crossing features 
of a light curve, compare Figure 9. We see from our results in 
Table 8 that the discovery of a moderately massive moon (see 
standard scenario assumptions in Table 1) is impossible for a 
source star size of /^source = 10 Rq or larger within our assump- 
tions. 

5.2.7. Increased separation between star and planet 

For reasons given above in Section 4.3, we did not vary the sep- 
aration between star and planet significantly, but stayed in the 
outer ranges of the planet lensing zone, corresponding to a well- 
detectable, wide-separation planetary caustic. For PS = 1.46 E , 
we get a detectability of 33.3 %, which is similar to our standard 
case (6 PS = 1.3 6 E ) with 30.6 %. 

6. Conclusion and outlook 

In this work we provide probabilities for the detection of moons 
around extrasolar planets with gravitational lensing. We showed 



that massive extrasolar moons can principally be detected and 
identified via the technique of Galactic microlensing. From our 
results it can be concluded that the detection of an extrasolar 
moon - under favourable conditions - is within close reach of 
available observing technology. The unambiguous characteri- 
sation of observed lunar features however will be challenging. 
The examined lens scenarios model realistic triple-lens config- 
urations. An observing rate of about one frame per 15 minutes 
is desirable, which is high, but well within the range of what is 
now regularly performed in microlensing follow-up observations 
of anomalous events. Similarly, an observational uncertainty of 
about 20 mmag can be met with today's ground-based telescopes 
and photometric reduction techniques for sufficiently bright tar- 
gets. Bright microlensing targets are mostly giant stars, which 
is an impediment to the detection of moons: Bulge giants with 
radii of order 10 R Q or larger smooth out any lunar caustic fea- 
ture. Therefore, in order to find moons, it is crucial to be able 
to perform precise photometry on small sources with angular 
sizes of the order of 10 -3 Einstein radii, corresponding to a few 
solar radii or smaller at a distance of 8 kpc. This means dwarf 
stars rather than giants need to be monitored in order to identify 
moons in the intervening planetary systems. Some improvement 
in resolution can be gained with the lucky-imaging technique on 
medium-sized ground telescopes (Grundahl et al. 2009). Under 
very favourable circumstances, exomoons might already be de- 
tectable from ground. However, future space-based telescopes, 
such as ESA's proposed Euclid 3 mission (cf. Refregier et al. 
(2010), Chapter 17) or the dedicated Microlensing Planet Finder 
(MPF) mission proposed to NASA (cf. Bennett et al. (2009)), 
will surely have the potential to find extrasolar moons through 
gravitational lensing. 
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ity density function /•(*,•) with a standard deviation of <j\. (cr]) 2 
is the variance of the distribution. 4 

We now want to examine whether the X t could be described 
equally well with a binary-lens model (J?. We introduce the new 
^-distributed 5 random variable 

1=1 V U i 

At this step, we could simulate data X t in order to find a some- 
what representative, randomly drawn value of Q 2 , but instead we 
simply calculate what the mean value of all possible Q 2 would 
be. We use the definitions, to find 



«*■>-£ 

> i=l 



Here, we keep in mind that, in general, (x 2 ) + (x) 2 . Using the 
parameters // and cr'. of the distribution /i(x,), we reduce the 
equation by calculating 

n 1 C°° 

i=l J-QO 
ft 1 r*oo 



Appendix A: Significance of deviation 

This section explains the method that we use to determine the 
significance of deviation between the two simulated light curves 
of the triple lens and the best-fit binary lens model. 

When we want to decide whether real observational data is 
better described with a triple-lens or with a binary-lens model, 
we can fit both models to the data and use the ^ 2 -test to assess 
the significance of the deviation between the data and the models 
and pick the model with the better fit. 

We do not have observational data, but numerically com- 
puted light curves. If we observed one of the simulated triple- 
lens systems, we would expect the data points to have a Gaussian 
distribution around the theoretical values. We could produce arti- 
ficial data by randomly scattering the simulated triple-lens light 
curve data points around their theoretical values and then fit a 
binary-lens model and the triple-lens model to this artificial data 
set, as described in Section 3.4. This approach allows us to use 
the usual ^ 2 -test, but is computationally intensive, because the 
models have to be fitted numerous times to get a statistically 
sound sample of ^-values. In our approach, we directly calcu- 
late the ^f 2 -value which is the expectation value of the above 
method. 

We hypothesise, every simulated point on our triple-lens 
light curve is in fact the mean // of the distribution of a random 
variable X t that is distributed according to the Gaussian probabil- 



We used the definitions of yuj and (crj) 2 and the property of the 
probability density function J fi(xi)dxi — 1 . In our case cr' = cr* 
holds true without loss of accuracy, and indeed, we simplify 
cr, = cr for all i, that is we assume that the photometric uncer- 
tainty is the same for all data points, and argue that this does not 
pose a problem as long as cr is chosen to match the maximum 
photometric uncertainty. So we can further reduce to 



<e 2 > 



n j ft 

i=l i=\ 



4-4 



This is the mean of all Q 2 possibly resulting, when compar- 
ing the simulated binary-lens light curve to randomly scattered 
triple-lens light curve points. It shall be our measure of devia- 
tion. Fortunately, we can provide all remaining parameters: 



4 We recall the definitions of mean and variance. For the random vari- 
able X with probability density function f(x), they are in general 

ji = I x f(x) dx = (X) and 

J— CO 

cr 2 = (x-fif f(x)dx = ((X-p)\ 

J —CO 

where we also introduce the notation (X) for the expectation value of X. 

5 The ^-distribution can be applied to the sum of n independent, 
normally distributed random variables Z ; with the mean flj = and the 
standard deviation <X; = 1 for all (', 

n 

See any standard book with an introduction to calculus of probability. 
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n is the number of degrees of freedom, equal to the number of 

compared data points, 
cr is the assumed standard deviation or photometric uncertainty 

of observations. 

(/ij - n 1 !) 2 equals the difference between the two compared light 
curves squared, which we already use for our least square fit. 



X 2 probability density function 





Ax 2 ) = degrees of freedom 




\<G 2 > 




^^PJ 0.01 





Fig.A.l. The^f 2 probability density function plotted for illus- 
tration only. For increasing degrees of freedom the central limit 
theorem takes effect and the curve will very closely resemble a 
Gaussian distribution. If P(X > (Q 2 )) < 1%, i.e. if the proba- 
bility P for a ^-distributed random variable X to be larger or 
as large as the mean of Q 2 is less than 0.01, we consider the 
deviation between the two compared curves to be significant. 



avoid refitting the binary-lens model light curve and always com- 
pare to the one we got as best-fit to the exact (simulated) triple- 
lens model light curve. This overestimation effect vanishes for 
n — > oo and is negligible for a sufficiently large n. The high sam- 
pling frequency we use, ensures that the number of data points 
is always large enough (n > 250 in all cases). 

While we cannot exactly quantify the uncertainty of the re- 
sults, we know that they pose strict upper limits for the lunar 
detectability in the various scenarios. This knowledge would en- 
able us to infer a tentative census of extrasolar moons, once the 
first successful microlensing detections have been achieved. 



Now we have to decide whether the deviation between our 
two examined light curves is significant. By taking (Q 2 ) as y 2 ^ 
one can evaluate the cumulative distribution function F„(z) = 
f-oo f( x )dx of the ^-distribution - we relied on the GNU 
Scientific Library 6 - to find the corresponding probability for any 
given^ 2 -distributed random variable X to be as large as or larger 
than <g 2 >, 

P(X>(Q 2 }) = l-FMQ 2 )). 

This probability is just the integral of the^f 2 probability density 
function to the right of (Q 2 ) as illustrated in Figure A.l. If this 
probability is very small, (Q 2 ) is outside the expected range for a 
X 1 -distributed random variable Q 2 . In that case, Q 2 is obviously 
not ^-distributed with n degrees of freedom, so we must con- 
clude that there is a significant deviation between the triple lens 
model light curve and the binary-lens model light curve. 

We say we have a significant deviation between our two 
curves, if the mean value (Q 2 ) is so high that the probability 
P for any random variable X being larger or equally large is 
less than 1%. We interpret this to say, only if the probability 
for a given triple lens light curve with independent, normally 
distributed data points to be random fluctuation of the compared 
binary-lens light curve is less than 1%, we consider it to be prin- 
cipally detectable. 

There are two known sources for overestimating the detec- 
tion rate coming with our method. First, systematic errors are 
not accounted for. A possible solution could be to add a further 

term to (Q 2 ) as in (Ql ew ) = {Q 2 )-n (^) , where cr sys can be as- 
sumed to lie in the few percent region for real data. Secondly, we 



6 http : //www . gnu . org/sof tware/gsl/ 



